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ABSTRACT 


An  annular  turbine  cascade  was  designed  and  manufactured  for  laser  Doppler  velocimetry 
and  probe  measurements  in  a  small-scale  rig.  The  purpose  of  the  experiment  was  to  determine 
the  limitations  of  these  measurements  in  a  confined  annulus,  and  to  compare  the  experimental 
results  with  numerical  predictions.  Downstream  probe  surveys  were  conducted  at  a  Reynolds 
number  of  500000.  Total  pressure  and  temperature  measurements  were  taken  upstream  of  the 
test  section,  and  hub  static  pressure  was  measured  downstream  of  the  cascade.  A  two- 
dimensional  laser  anemometer  was  used  to  obtain  preliminary  velocity  field  measurements  for 
code  validation.  The  two  main  problem  areas  identified  were  the  type  of  seeding  material  used 
and  the  configuration  of  the  optical  access  window.  Flowfield  characteristics  were  predicted 
using  a  three-dimensional  viscous  code,  which  were  then  compared  to  experimental 
measurements.  The  comparison  showed  that  the  numerical  simulator  predicted  well  the  general 
features  of  the  flow  field  that  were  measured  in  the  experiment.  Recommendations  are  made 
which  would  improve  the  mapping  of  the  velocity  field  for  a  more  complete  code  validation. 
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I .  INTRODUCTION 


The  measurement  of  secondary  flows  in  annular  turbine 
cascades  (ATC)  is  very  important  to  gas  turbine  designers 
and  computational  fluid  dynamicists.  A  decrease  in  stator 
performance  caused  by  secondary  flow  has  a  significant 
effect  on  turbine  stage  performance.  Therefore,  the 
quantification  of  secondary  flow  effects  on  performance 
(i.e.  losses)  is  particularly  important  for  small  aspect 
ratio  core  turbines. 

The  computational  prediction  of  component  performance  is 
potentially  a  valuable  tool  for  turbomachinery  designers.  A 
small  improvement  in  engine  efficiency  can  amount  to  huge 
savings  in  yearly  fuel  costs  for  a  wing  of  aircraft.  Good 
quantitative  predictions  of  turbomachinery  performance  are 
necessary  for  turbomachinery  designers  to  increase  engine 
performance  at  a  minimum  cost. 

This  report  presents  the  results  of  numerical  simulation 
of  an  ATC  and  verification  by  experimental  measurements.  A 
cascade  was  designed  and  built  by  the  author  so  that  laser 
Doppler  velocimetry  (LDV)  and  probe  measurements  of  turbine 
secondary  flows  could  be  performed.  Three- hole  probe 
measurements  were  taken  in  the  wake  of  a  high  turning 
turbine  nozzle.  These  were  compared  with  the  numerical 


results  obtained  by  running  RVC3D,  a  three-dimensional, 
viscous,  computer  program. 

The  longer  term  goal  of  this  research  is  to  obtain  LDV 
measurements  as  close  as  possible  to  the  endwalls  and  to 
fully  determine  the  flowfield  structure  in  the  presence  of 
tip  clearance  effects. 

Several  researchers  have  investigated  secondary  flows  in 
turbine  cascades.  Louis  Goldman  and  Richard  G.  Seaholtz 
have  published  several  papers  on  laser  measurements  in  an 
annular  cascade  of  high  turning  core  turbine  vanes  [Ref.  1  & 
2] .  Their  work  consisted  of  making  laser  anemometer 
measurements  in  the  blade  passages,  and  limited  measurements 
in  the  wake  region.  These  measurements,  along  with  blade 
surface  static  pressures,  were  compared  to  predictions  made 
with  a  three-dimensional  inviscid  flow  analysis  prograun. 

The  experimental  measurements  compared  well  with  the 
calculated  values.  There  was  poor  agreement  however  in  the 
wake  region  where  viscous  effects  dominated,  which  could  not 
be  resolved  with  the  inviscid  flow  analysis  program. 

A.  Yamamoto  ha.«)  also  published  several  papers  on 
secondary  flows  in  turbine  cascades.  Reference  3  reports  a 
detailed  investigation  of  secondairy  flow/loss  mechanisms. 
Five -hole  probe  measurements  in  two  types  of  turbine 
cascades  with  different  turning  angles  were  presented. 
Ycuncunoto's  work  focused  on  complete  flowfield  surveys  within 
the  blade  rows  of  linear  turbine  cascades. 
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Rodrick  V.  Chima  developed  a  three-dimensional  flow 
analysis  code  to  compute  the  design,  and  off -design 
operating  conditions  in  transonic  rotors  [Ref.  4).  His 
code,  RVC3D,  a  computer  code  for  analysis  of  three- 
dimensional  viscous  flows  in  turbomachinery,  was  used  for 
the  numerical  prediction  of  the  flows  in  the  annular  turbine 
cascade . 
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II.  EXPERIMENTAL  SETUP 


A.  APPARATUS 

The  experimental  results  were  obtained  using  a  transonic 
turbine  stator  tested  in  an  annular  cascade  facility  at  the 
Turbopropulsion  Laboratory  of  the  Naval  Postgraduate  School. 
The  schematic  of  the  test  rig  can  be  seen  in  Figure  1.  The 
hub  and  tip  have  constant  radii  with  a  hub  radius  of  3.89  in 
and  a  tip  radius  of  4.585  in.  The  tested  cascade  had  31 
blades  measuring  1.1  in  high  and  1.00  in  axial  chord,  as  can 
be  seen  in  Figure  2.  The  blade  spacing  at  mid -span  was  0.86 
inches.  The  blades  were  stacked  radially,  and  there  was  no 
change  in  blade  shape  in  the  radial  direction.  The  inflow 
was  purely  axial  at  an  average  Mach  number  of  0.13.  The 
inlet  total  -  to  -  exit  hub  static  pressure  ratio  was  0.68, 
which  gave  an  average  exit  Mach  number  of  0.59.  The 
apparatus  that  housed  the  turbine  can  be  seen  in  Figures  3 
and  4  . 

Airflow  was  provided  by  a  VA-312  Allis  -  Chalmers ,  12- 
stage  axial -flow  compressor,  operating  at  12,000  rpm.  The 
compressor  produced  a  mass  flow  rate  of  7.74  Ib/sec,  during 
the  experiment.  The  piping  schematic  can  be  seen  in  Figure 
5.  The  air  entered  the  ten  inch  flange  (Fig.  4)  through  a 
honeycomb  flow  straightener ,  after  leaving  the  plenum 


INLET 

FLOW 


Plgur«  2.  Blade  geometry 
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chamber.  The  air  was  then  directed  into  the  test  section 

through  the  bellmouth,  (Fig.  l)  which  smoothly  reduced  the  • 

airflow  from  a  ten  inch  diameter  to  a  9.17  inch  diameter 

section.  After  leaving  the  turbine  blades  the  air  was 

dvimped  three  chord  lengths  beyond  the  blade  row,  on  the  • 

inner-hub  surface.  The  air  was  finally  discharged  at  the 

outer  tip  surface,  which  was  four  chord  lengths  downstream, 

into  the  test  cell.  • 

Laser  measurements  were  designed  to  be  taken  through  an 
optical  window  which  was  centered  two  chords  downstream  from 
the  trailing  edge.  As  seen  in  Figure  6,  the  LDV  measurement  • 

points  were  located  downstream  of  the  blades  at  different 
peripheral  stations.  These  tangential  surveys,  at  different 
radial  locations,  would  cover  a  plane  of  data  output  from  • 

the  numerical  simulation.  This  was  to  allow  direct 
comparison  between  the  numerical  simulation  and  the 

experimental  results.  • 

B.  COBRA  PROBE  MEASUREMENTS 

A  cobra  probe  was  used  to  measure  the  flowfield  on  the 
exit  plane  at  various  radii,  and  at  the  same  axial  location 
as  the  LDV  measurements.  The  cobra  probe  measurements  were 
taken  at  10,  25,  50,  75,  and  100  percent  span,  as  measured 
from  the  hub. 

The  cobra  probe  was  a  three-hole  probe,  with  the  center 
hole  0.036  Inches  in  diameter,  and  the  two  side  ports  set  at 


•  • 


•  •  •  •  • 


10 


an  angle  of  50  degrees.  The  probe  was  calibrated  in  a  free 
jet  in  which  a  pitot  static  probe  was  used  to  obtain  the 
Mach  number  of  the  flow.  This  was  plotted  against  the 
indicated  total  head,  to  obtain  a  correlation  for  the  cobra 
probe.  The  results  can  be  seen  in  Appendix  C. 

C.  LASER  SETUP 

A  TSI  Inc.  four-beam,  two  color  fiber-optic  LDV  system 
was  set  up  and  used  as  a  one-component  system  in  backscatter 
mode.  Components  of  the  LDV  optics  are  shown  in  Figure  7, 

The  laser  was  a  Lexel  Model  95  four-Watt  argon-ion  laser 
operating  nominally  at  two  Watts.  The  LDV  system  used  four 
argon-ion  laser  beams,  two  blue  (488  nm)  and  two  green 
(514.5  nm)  from  a  modular  color  separator,  to  measure  two 
orthogonal  velocity  vectors.  The  color  separator,  and  its 
frequency  shifting  Bragg  cell,  comprised  the  "Colorburst" . 
The  beams  were  coupled  into  four  single-mode  polarization¬ 
preserving  fibers  using  two  translator  modules  and  two 
double-input  couplers. 

The  optical  fibers  carried  the  laser  power  to  the  probe 
head,  where  the  light  from  each  fiber  was  collimated.  The 
four  collimated  beams  were  focussed  and  crossed  through  a 
transmitting  lens.  Backscattered  blue  and  green  light  was 
collected  through  the  transmitting  and  receiving  lenses  and 
focused  into  the  multimode  receiving  fiber.  The  received 
light  from  the  multimode  fiber  was  collimated  with  the 


Figure  7 
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double- input  coupler.  The  collimated  received  light  was 
directed  into  a  color  separator  to  separate  the  light  from 
the  two  velocity  components.  This  color  separator  and  its 
photomultipliers  and  frequency  shifters  comprised  the 
"Col or link" . 

LDV  signals  were  processed  by  two  TSI  Model  1990C 
counter- type  signal  processors  as  shown  schematically  in 
Figure  8.  The  Model  1990C  measured  the  time  taken  for  a 
given  number  of  cycles  in  a  Doppler  burst  using  a  high 
resolution  clock  (±1  ns) . 

An  oscilloscope  attached  to  the  input  conditioner 
provided  real-time  display  of  the  photomultiplier  output  for 
setting  filters  on  the  counters.  The  counters  were  operated 
in  a  single  measurement  per  burst  (SM/B) ,  coincident  mode. 

A  digital  interface  on  the  counter  provided  two  functions. 
First,  the  master  interface  compared  the  incoming  signal 
from  each  counter  and  checked  for  coincidence  validation. 
Second,  the  interface  provided  computer  input  using  direct 
memory  access,  sending  five  16  bit  words  for  each  valid 
burst  to  the  computer. 

The  full  three-dimensional  set-up  and  description  of 
this  system  is  given  in  Reference  5. 
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TEST  SECTION 


D.  EXPERIMENTAL  PROCEDURE 

When  supply  conditions  were  steady,  measurements  were 
taken  during  the  e'cperiment  that  were  used  as  input  for 
RVC3D.  These  included;  upstream  temperature  and  pressure, 
and  downstream  static  pressure  at  the  hub.  These  parameters 
were  used  as  inputs  to  the  computer  program  for  the 
numerical  simulation  of  the  flow  through  the  ATC. 

Flowfield  measurements  were  taken  two  axial  chord 
lengths  downstream  from  the  trailing  edge  of  the  vanes,  with 
a  cobra  probe.  A  ruler  was  used  to  position  the  probe  when 
it  was  rotated  to  the  next  peripheral  measurement  location. 

Initially,  course  surveys  were  performed  across  one 
blade  passage.  Next,  finer  circumferential  surveys  were 
performed  across  two  blade  wakes.  This  gave  better 
resolution  of  the  wake  profiles,  and  also  a  check  for  flow 
periodicity.  The  initial  coarse  surveys,  at  0.1  inch 
spacing  were  completed  at  25,  50,  and  75  percent  span.  The 
graphical  representation  of  these  data  can  be  seen  in 
Appendix  C.  The  finer  set  of  measurements  was  taken  at  10, 
25,  50,  75,  and  90  percent  span,  spaced  at  0.05  inches.  The 
measurements  were  also  extended  to  capture  two  wake  regions, 
to  check  for  periodicity.  The  tabulated  data  are  given  in 
Appendix  C. 
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III.  NUMERICAL  SIMULATION 


A.  COMPUTATIONAL  SCHEME 

RVC3D  (Rotor  Viscous  Code  3-D)  is  a  computer  code  for 
the  analysis  of  three-dimensional  viscous  flows  in 
turbomachinery.  The  code  solves  the  thin- layer  Navier- 
Stokes  equations  with  an  explicit  finite-difference 
technique.  It  is  applicable  to  annular  blade  rows  or  linear 
cascades.  Two  algebraic  turbulence  models  and  a  simple  tip 
clearance  model  are  available. 

The  code  solves  the  Navier- Stokes  equations  formulated 
in  a  Cartesian  coordinate  system  with  rotation  about  the  x- 
axis.  The  equations  are  mapped  to  a  general  body- fitted 
coordinate  system.  Streamwise  viscous  terms  are  neglected 
using  the  thin- layer  assumption,  but  all  cross  -  channel 
viscous  terms  are  retained.  Turbulence  effects  are  modeled 
using  either  a  3-D  adaptation  of  the  Baldwin- Lomax 
turbulence  model  [Ref.  6]  or  the  Cebeci -Smith  model  [Ref. 

7] .  The  equations  are  discretized  using  second-order 
finite-differences  and  solved  using  a  multistage  Runge-Kutta 
scheme.  References  4  and  8  describe  the  mathematical 
formulation  of  the  RVC3D  code.  A  saunple  input  file  for 
RVC3D  is  presented  in  Appendix  A. 


•  • 
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B .  GRID  GENERATION 


The  grid  generation  started  with  the  schematic  the 
ATC,  as  seen  in  Figure  1.  The  drawings  were  digitized,  to 
produce  a  set  of  x-y  points.  These  data  points  were  used  as 
part  of  the  input  file  for  the  GRAPE  code  iRef .  9] .  The 
input  file  for  the  GRAPE  code  is  presented  in  Appendix  A. 
This  code  produced  a  two-dimensional  C-type  grid  about  the 
turbine  blade,  as  seen  in  Figure  9.  Another  FORTRAN  code, 
STACK,  takes  a  two-dimensional  grid  and  stacks  it  into 
three  -  dimensions ,  using  hyperbolic  tangent  stretching. 

STACK  produces  a  grid  output  file  that  can  be  used  by  RVC3D. 
The  namelist  input  for  STACK  is  given  in  Appendix  A.  The 
three-dimensional  grid,  which  has  121x31x21  grid  points,  can 
be  seen  in  Figure  10. 

A  second  grid  generation  program  was  also  utilized. 
TCGRID  {Turbomachinery  C  GRID)  generates  three-dimensional 
C-  or  H-type  grids  for  turbomachinery.  TCGRID  generates  a 
computational  grid  in  a  single  step  once  the  blade  surfaces 
have  been  defined.  The  x-y  coordinates,  from  the 
digitization  of  the  drawings,  were  translated  to  z-r- theta 
coordinates  that  TCGRID  could  use.  A  FORTRAN  program  was 
written  to  accomplish  this,  and  it  is  included  in  Appendix 
A.  Once  the  coordinates  were  in  the  proper  format  they  were 
included  in  the  input  file  for  TCGR'^0,  which  then  produced 
an  output  file  compatible  with  RVC3D.  The  input  file  for 
TCGRID  is  presented  in  Appendix  A. 


Figure  10.  3-D  grid 


TCGRID  generates  grids  using  the  following  technique: 

1.  A  coarse,  equally- spaced  meridional  grid  is  generated 
between  the  supplied  hub  and  tip. 

2.  The  blade  geometry  is  interpolated  onto  the 
meridional  grid. 

3.  2-D  blade- to -blade  grids  are  generated  along  the 
meridional  grid  lines  in  (m,  rbar*theta)  coordinates, 
using  a  version  of  the  Sorenson  [9]  GRAPE  code. 

4.  The  (m, rbar*theta)  coordinates  are  transformed  to 
( z , r , theta) . 

5.  The  2-D  grids  are  reclustered  spanwise  to  make  a  full 
3-D  grid. 

6.  Finally  the  (2,r, theta)  coordinates  are  transformed 
back  to  (x,y,z)  and  written  in  PL0T3D  format. 

When  the  final  grid  was  complete  it  was  stored  in  a 
format  that  RVC3D  could  use  to  run  the  flowfield  solution. 
After  the  solution  was  obtained,  it  was  analyzed  to  ensure 
that  it  was  grid  independent.  One  way  to  check  for  grid 
independence  is  that  the  wake  should  follow  the  trailing 
edge  projection  and  not  the  grid  centerline. 


* 


4r 


1.  GRID  GENERATION  COMPARISON 

The  grid  generation  programs  both  produced  an  * 

acceptable  grid.  The  main  difference  between  the  two  was 
that  TCGRID  was  a  one -step  process,  while  STACK  required  a 
GRAPE  grid  as  an  input  file.  It  took  two  steps  to  produce  a  * 

final  grid  using  STACK,  while  TCGRID  could  produce  it  in  one 
step.  TCGRID  also  had  more  options  available  to  fine-tune 
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the  grid.  GRAPE/STACK  was  initially  easier  to  use,  due  to 
prior  experience  with  GRAPE,  which  allowed  the  first  grids 
to  be  produced  quickly. 

TCGRID  required  (z.r, theta)  coordinates,  which  were 
not  readily  available  from  the  drawings.  After  data 
conversion,  TCGRID  proved  very  useful,  and  the  easier  of  the 
two  to  manipulate. 

2 .  RBCOMMEMDATIOKS 

If  the  turbine  geometry  data  are  in  x-y  format,  the 
GRAPE/STACK  combination  is  quite  effective  for  producing 
early  grids.  These  grids  can  be  used  to  test  the  flowfield 
solver.  GRAPE/STACK  grid  generation  occurs  in  discrete 
steps,  which  helps  in  the  process  of  fine-tuning  the  final 
grid.  One  drawback  is  that  STACK  will  only  stack  radially. 

TCGRID  should  be  used  if  the  turbine  coordinates  are 
given  in  a  format  that  TCGRID  can  accept.  TCGRID  would  also 
be  recommended  since  it  is  a  one  step  grid  generation 
progreim  that  has  powerful  tools  for  manipulating  the  final 
grid.  TCGRID  can  account  for  twist  of  the  airfoil,  as  well 
as  non- radial  stacking  and  end -bends. 


IV.  RESULTS  AND  DISCUSSION 


A.  LASER  DOPPLER  VELOCIMETRY  MEASUREMENTS 

Three  different  materials  were  use  for  seeding;  olive 
oil,  water,  and  a  glycerine  oil/water  mixture.  The  olive 
oil  gave  the  best  LDV  signal,  but  contaminated  the  window  at 
higher  speeds.  The  final  LDV  configuration  was  established 
by  directing  the  probe  into  the  wake  region  from  downstream 
of  the  cascade,  as  can  be  seen  in  Figure  11.  Because  of  the 
current  geometry  of  the  test  section  only  the 
circumferential  velocity  was  measured.  A  histogram  of  the 
LDV  output  is  given  in  Appendix  F.  The  measured  flow  of 
this  histogram  was  76.695  M/S,  at  a  turbulence  intensity  of 
6.488%.  At  speeds  higher  than  this,  seeding  became  the 
primary  problem. 

B.  NUMERICAL  RESULTS 

Figure  12  shows  the  computed  blade  surface  static 
pressure.  The  figure  shows  results  that  were  expected.  The 
pressure  is  higher  on  the  tip  than  at  a  corresponding  point 
on  the  hub.  The  suction  side  has  a  steeper  pressure 
gradient  than  the  pressure  side,  and  the  pressures  are  equal 
at  the  trailing  edge.  Figure  13  shows  the  normalized  blade 
static  pressure  at  mid-span.  The  discontinuity  on  the 
suction  side  is  due  to  the  sonic  region.  Figure  14  shows 


PRESSl  RF 


Blade  surface  static  pressure 


Figure  13.  Blade  surface  static  pressure  at  mid- span 
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Figure  14.  Stagnation  pressure  on  the  exit  plane 


the  stagnation  pressure  on  a  exit  plane  one  chord  length 
from  the  trailing  edge.  The  wakes  are  clearly  evident  in 
this  figure.  Unfortunately  the  computational  grid  could  not 
be  extended  two  chord  lengths  downstreaun  of  the  trailing 
edge.  Attempts  at  producing  larger  grids  gave  rise  to 
unacceptable  shapes.  All  the  flow  field  comparisons  with 
the  experiment  were  made  using  data  computed  at  one  chord, 
whereas  the  experimental  data  were  measured  at  two  chord 
lengths  downstream. 

C.  COMPARISON  OF  DOWNSTREAM  LOSSES  AMD  EXIT  FLOW  ANGLES 

Figure  15  is  a  comparison  of  the  Mach  ntunber  at  10% 
span  predicted  by  the  computer  code,  with  the  experimental 
results  from  the  cobra  probe.  The  velocity  peak,  at  the  20% 
chord  point,  appears  to  be  a  jet,  but  upon  further  analysis 
it  appears  that  the  hub  boundary  layer  is  being  pinched,  and 
the  freestream  core  flow  is  migrating  toward  the  hub 
surface.  This  can  be  clearly  seen  in  Figure  14  where  the 
blue  region  is  thinnest  close  to  the  hub  (i.e.  the  core  flow 
is  migrating  toward  the  hub  endwall) . 

Figures  16,  17,  18,  and  19  show  the  Mach  nximber 
comparisons  for  flow  at  25,  50,  75,  and  90  percent  span 
respectfully.  The  computer  code  does  a  good  job  predicting 
the  wake  deficit,  and  the  width  of  the  wake.  There  is  good 
agreement  at  all  radial  locations.  The  core  flow  Mach 
number  is  about  0.65  and  the  jet  flow  in  Figure  15  is  only 
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Figure  15 .  Mach  number  coiig)ari8on  at  10%  span 
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Figure  16.  Mach  number  comparison  at  25%  span 
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Figure  18.  Mach  ntunber  comparisons  at  ISH  span 


•  • 


Figure  19.  Mach  number  comparisons  at  90%  span 


I 


at  a  Mach  number  of  0.5.  This  suggests  a  rapid  growth  of 
the  hub  endwall  boundary  layer  rather  than  a  jet  flow. 


Figures  20-24  shows  the  flow  angle  comparisons  for  the 
five  different  radial  positions  that  were  surveyed.  The 
overall  trend  in  flow  angle  variation  in  the  circumferential 
direction  is  predicted,  but  the  magnitude  of  the  flow  angle 
is  not.  The  computational  grid  may  need  to  be  reconfigured 
with  an  exit  angle  closer  to  the  exit  flow  angle  in  an 
attempt  to  improve  the  comparison,  once  the  axial  locations 
of  experiment  and  computation  coincide. 
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Figure  20.  Flow  angle  comparison  at  10%  span 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

Probe  measurements  were  easy  to  take  and  were  reliable. 
They  allowed  a  quick,  qualitative  evaluation  to  be  made  of 
the  computer  code.  The  spatial  resolution  and  accuracy  were 
not  as  good  as  should  be  expected  from  the  laser  system. 

A  further  problem  with  the  cobra  probe  was  in  the 
measurement  of  the  flow  angle.  A  hysteresis  effect  was 
evident  depending  on  which  way  the  probe  was  turned.  This 
was  compensated  for  by  overshooting  the  new  angle  and  then 
returning  to  the  pneumatic  balance  point.  The  settling  time 
on  the  yaw  manometer  was  very  slow  due  to  the  size  of  the 
static  ports. 

The  computer  model  qualitatively  predicted  the  measure 
of  variation  in  the  mean  flow  Mach  number.  The  flow  angle 
variation  prediction  was  reasonable,  given  that  it  is  a 
second  order  effect. 

The  computer  model  had  run  times  of  four  hours  on  the 
Stardent,  and  12  hours  on  the  IRIS  workstations.  Clearly, 
all  calculations  should  be  done  locally  on  the  Stardent,  and 
when  resources  are  made  available,  on  the  Cray  YMP  at  the 
National  Aerodynamic  Simulator. 

In  order  to  get  full  three-dimensional  laser  data,  a 
seeding  material  needs  to  be  found  that  will  not  dirty  the 
optical  access  window,  but  will  follow  the  streamlines. 


Water  evaporated  before  reaching  the  probe  volume,  and 
solids  were  not  used  due  to  their  possible  toxic  nature  in 
an  enclosed  facility.  Olive  oil  was  used,  but  the  optical 
access  window  became  dirty  almost  instantly. 

The  second  hurdle  that  needs  to  be  overcome  is  the  glare 
caused  by  the  laser  beams  reflecting  off  the  hub.  The 
reflections  quickly  saturate  the  photomultiplier,  greatly 
reducing  the  signal  to  noise  ratio.  Painting  the  hub  black 
helped  reduce  the  noise  significantly.  Another  reduction 
technique  was  to  put  an  opaque  screen  between  the  probe  and 
the  probe  volume.  The  screen  had  cutouts  to  allow  the  laser 
beams  and  the  return  signal  clear  passage,  but  blocked  all 
other  reflections. 

A  possible  solution,  to  the  glare  problem,  is  to  have 
the  laser  entering  the  window  at  a  high  angle  so  that  the 
beams  are  reflected  away  from  the  optical  window.  This  was 
attempted  on  a  limited  basis,  but  a  more  thorough  effort 
needs  to  be  undertaken.  A  second  alternative  is  to  attempt 
to  get  laser  measurements  by  going  back  into  the  test  cell 
from  a  downstream  location.  This  would  generate 
significantly  less  glare.  The  last  ring  of  the  apparatus 
will  need  to  be  modified  to  allow  full  optical  access.  Once 
this  modification  is  made,  measurements  could  be  taken  right 
up  to  the  trailing  edge.  With  these  changes  LDV  data  could 
be  taken  close  to  the  endwall,  and  the  flowfield  determined 
in  the  presence  of  tip  clearance  effects. 


« 
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APPENDIX  A.  SAMPLE  INPUT  PILES 


This  file  is  the  input  file  for  TCGRID,  for  the  annular 
cascade  turbine. 

&naml  iin=121  jm=31  kni=21  itl=21  icap=12  k2d=3  merid=0  &end 
&nani2  nle=10  nte=lO  dsle=.020  dste=.003  dshub=.  00004 
dstip=. 00004  dswte=.0003  dswex=.040  dsthr=l. 

dsinin=.0004  dsmax=.025 

dsra=.45  rcorn=.098  &end 

&nam3  itenn=l00  idbg=0  0000000  aabb*.5  ccdd=.45  &end 
&nain4  zbc=-1.5  -1.5  1.0  -1.5  -1.5  1.0 

rbc=  3.895  3.895  3.895  4.585  4.585  4.585  &end 
'TRANSONIC  TURBINE' 

2  2 

-.975  0.0 

3.895  3.895 

-.975  0.0 

4.585  4.585 

2  51  31 

-7. 8999996E-03  - 9 . 9200001E- 03  - 1 . 2000000E- 02  - 1 . 4080000E- 02 
-1.6100001E-02  -1.7999999E-02  - 1 . 9710001E- 02  - 2 . 1190001E- 02 
-2.2390001E-02  -2 . 3280000E-02  -  8 . 1040002E- 02  -0.1503800 
-0.2308900  -0.3219400  -0.4228500  -0.5328500  -0.6511000 

-0.7893000  -0.9275000  -0.9410600  -0.9527700  -0.9622700 

-0.9692700  -0.9735600  -0.9750000  -0.9735600  -0.9692700 

-0.9622700  -0.9527700  -0.9410700  -0.9275000  -0.9124900 

-0.8250000  -0.7200000  -0.6000000  -0.4560000  -0.3440000 

-0.2730000  -0.2200000  -0.1820000  -0.1160000  -  5 . 9000000E- 02 

-1.7000001E-02  O.OOOOOOOE+00  - 1 . 8220000E- 04  - 7 . 2359998E- 04 
-1.6100000E-03  -2.8100000E-03  -4 . 2900001E- 03  -  6 . OOOOOOIE- 03 
-7. 8999996E-03  -0.2349885  -0.2351271  -0.2351733  -0.2351271 
-0.2349885  -0.2347600  -0.2344519  -0.2340719  -0.2336329 

-0.2331451  -0.1980565  -0.1643980  -0.1324210  -0.1023671 

-7.4464694E-02  -4 . 8926830E-02  -2 . 5946086E- 02  - 4 . 6726577E- 03 
1.6618744E-02  1 . 9057767E- 02  2 . 2061616E- 02  2 . 5545571E- 02 
2.9399229E-02  3 . 3504494E- 02  3 . 7740692E- 02  4 . 1976891E- 02 
4.6082158E-02  4 . 9935814E- 02  5 . 3417202E- 02  5 . 6423619E- 02 
5.8862645E-02  6 . 0659818E- 02  6 . 9319643E- 02  7 . 1887039E- 02 
6.7522466E-02  5 . 1347882E- 02  2 . 5673941E-02  0 . OOOOOOOE+00 
-2.5673941E-02  -  5 . 1347882E- 02  -0.1026958  -0.1540437 
-0.2053915  -0.2320924  -0.2326265  -0.2331451  -0.2336329 

-0.2340719  -0.2344519  -0.2347600  -0.2349885  3.895000 

3.895000  3.895000  3.895000  3.895000  3.895000  3.895000 

3.895000  3.895000  3.895000  3.895000  3.895000  3.895000 

3.895000  3.895000  3.895000  3.895000  3.895000  3.895000 

3.895000  3.895000  3.895000  3.895000  3.895000  3.895000 

3.895000  3.895000  3.895000  3.895000  3.895000  3.895000 

3.895000  3.895000  3.895000  3.895000  3.895000  3.895000 


3  . 
3  . 

3  . 
-1 
-1 
-8 
-0 

0 

-0 
-0 
-0 
-0 
-1 
-4 
-0 
-0 
-0 
-4 
1 . 
2  . 

4 

5  . 
4. 
-4 
-0 
-0 
4  . 
4. 
4  . 
4  . 
4  . 
4  . 
4. 
4  . 


095000  3. 
895000  3. 
895000  3. 
.2000000E 
.9710001E 
. 1040002E 
.5328500 
.9527700 
.9735600 
.9275000 
.4560000 
. 1160000 
.8220000E 
.2900001E 
. 1997426 
. 1991690 
. 1396576 
. 1563794E 
6189748E- 
8462378E- 
2420939E- 
1531076E- 
3620501E- 
.3620501E 
. 1971647 
.1991690 
585000  4. 
585000  4. 
585000  4. 
505000  4. 
585000  4. 
505000  4. 
585000  4. 
585000  4. 


895000 
-02 
-02 


095000  3.095000  3.895000  3.895000  3.895000 
895000  3.895000  3.095000  3.895000  3.895000 
-7.8999996E-03  -  9 . 9200001E- 03 
-1.4080000E-02  - 1 . 6100001E- 02  - 1 . 7999999E - 02 
-2 . 1190001E-03  -2 .2390001E-02  -  2 . 3280000E- 02 
-02  -0.1503800  -0.2308900  -0.3219400  -0.4228500 
-0.6511000  -0.7893000  -0.9275000  -0.9410600 
-0.9622700  -0.9692700  -0.9735600  -0.9750000 
-0.9692700  -0.9622700  -0.9527700  -0.9410700 
-0.9124900  -0.8250000  -0.7200000  -0.6000000 
-0.3440000  -0.2730000  -0.2200000  -0.1820000 
-5.9000000E-02  - 1 . 7000001E - 02  0 . OOOOOOOE  +  00 
-04  -7.2359998E-04  - 1 . 6100000E- 03  -  2 . 8100000E- 03 
-03  -6.0000001E-03  -  7 . 8999996E- 03  -0.1996249 
-0.1997019  -0.1997426  -0.1996249  -0.1994308 
-0.1988462  -0.1984733  -0.1980589  -0.1682508 
-0.1124929  -8.6961828E-02  -  6 . 3258447E- 02 
-02  -2 .2041440E-02  -  3 . 9694658E - 03  1 . 4 117776E - 02 
02  1 . 8741548E-02  2 . 1701200E- 02  2 . 4974918E- 02 
02  3 .2061070E-02  3 . 5659760E- 02  3 . 9147221E- 02 
02  4.5378406E-02  4 . 7932386E- 02  5 . 0004359E- 02 
02  5 . 8887679E-02  6 . 1068702E- 02  5 . 7360962E- 02 
02  2.1810250E-02  0 . OOOOOOOE+00  -2 . 1810250E- 02 
02  -8 . 7241001E-02  -0.1308615  -0.1744820 
-0.1976183  -0.1980589  -0.1984733  -0.1988462 
-0.1994308  -0.1996249  4.585000  4.585000  4.585000 
585000  4.585000  4.585000  4.585000  4.585000 


585000  4.585000  4.505000 
585000  4.585000  4.585000 
585000  4.585000  4.585000 
505000  4.505000  4.585000 
505000  4.585000  4.585000 
505000  4.585000  4.585000 


585000  4.585000  4.585000 


,585000 

,585000 

,585000 

,585000 

,585000 

.585000 

585000 


585000 

585000 

585000 

585000 

585000 

505000 

585000 


This  is  the  input  for  RVC3D: 


'TRANSONIC  TURBINE  Annular  Cascade' 

&nll  iin=l21  jm=31  kin=21  itl=21  iil=54  &end 

&nl2  cfl=5.5  aviscl=0.0  avisc2=0.0  avisc4=0.30  ivdt=l 

nstg=4  itmax=ll00  irs=l  epi=0.50  epj=0.60  ep)c=0.60 
&end 

&nl3  ibcin=3  ibcex=3  isymt=0  ires=10  icrnt=50 
iresti=0  iresto=l  ibcpw=0  iqin=0  &end 
&nl4  emxx=0.13  emty=0.0  emrz=0.0  expt=0.0  prat=0.6800 
ga=1.4  om=0. 000000  igeom=l  alex=-67.0  &end 
5cnl5  ilt=3  tw=l.00  renr=5.000e5  prnr=.7  prtr=.9 

vispwr= .  666666  srtip*0.0  cinutin=14.  jedge=15  kedge=ll 
iltin=2  dblh=0.0024  dblt=0.0089  &end 
&nl6  iol=l  io2=165  oar=0.  ixjb=0  njo=l  nko=3 
jo=l  ]co=5  11  16  send 
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This  is  the  input  file  to  GRAPE,  for  the  Annular  Turbine 
Cascade  : 

&GRID1  jmax=12l , kmax=21 , ntetyp=3 , nairf =5 , jairf =51 , nibdst=7, 
j  tebot=40 , j  tetop=121 , nobshp=7 , xleft  =  -1.5,xle=- .975, 
xte=0 . 0 , xright= . 75 , rcorn= . 08 , dsi= . 001 , norda=4 , 3 , 
maxita=200 , 400 , nout=4  &end 
&GRID2  dsra= . 45 , pitch= . 7895  &end 
&GRID3 

airfx=- .0079000, - .0099200, - .0120000, - .0140800, - .0161000, 

- .0180000, - .0197100, - .0211900, - .0223900, - .0232800, - .0810400, 
- .1503800, - .2308900, - .3219400, - .4228500, - .5328500, - .6511000, 
- .7893000, - .9275000, - .9410600, - .9527700, - .9622700,  -  .9692700, 
- .9735600, - .9750000, - .9735600, - .9692700, - .9622700, - .9527700, 
- .9410700, - .9275000, - .9124900, - .8250000, - .7200000, - .6000000, 
- .4560000,  -  .3440000,  - .2730000,  - .2200000, - .1820000,  - .1160000, 
- .0590000,  - .0170000,0.0000000, - .0001822, - . 0007236 ,  -  . 0016100 , 
-  .0028100,  - .0042900, - .00600000, - .0079000, 

airfy=- . 9152800 , - . 9158200 , - .9160000, - .9158200, - .9152800, 

- .9143900, -.9131900, -.9117100,- .9100000, - .9081000, - .7714300, 
- .6403300,  - .5157800, - .3987200,  - .2900400, - .1905700,  -  .1010600, 
- .0182000,0.0647300,0.0742300,0.0859300,0.0995000,0.1145100, 
0.1305000,0.1470000,0.1635000,0.1794900,0.1945000,0.2080600, 
0.2197700,0.2292700,0.2362700,0.2700000,0.2800000,0.2630000, 
0.2000000,0.1000000,0.0000000, - .1000000, - .2000000,  -  .4000000, 
- .6000000, - .8000000, - .9040000, - .9060800, - .9081000, - .9100000, 
- .9117100, - .9131900, - .9143900, - .9152800, 

&end 


This  is  the  namelist  for  STACK: 

&nll  km=21  rhub=3.895  rtip=4.585  nblade=31  ysp=-.3  dhl=0.01 
dtl=0.01  &end 


k=«lck(3) 
i  =  l 

do  10  j  =jmax, 1,-1 
um=uin+l 
do  15  1=1,5 
k=kk(l) 

rho=qq ( 1 , i , j  ,  k) 

v(l) = ( (qq{2. i, j ,k) /rho) **2+ (qq(3 , i. j , k) /rho) **2) **.5 
ang (1) =atan (qq (3 , i, j , k) /qq (2 , i,j,k))*57.3 
15  continue 
uinper=um/62 

write (7, 300)umper, (v(l) ,1=1,5) 
write (8, 300) umper, (ang(l) ,1=1,5) 

10  continue 
i  =  121 

do  30  j=l,jmax,l 
um=um+l 
do  25  1=1,5 
k=kk(l) 

rho=qq(l,  i,  j  ,)c) 

v(l)  =  (  (qq(2,i,  j,k)/rho)**2+(qq(3,i,  j,)c)/rho)**2)**.5 

ang  (1)  =atan  (qq  (3  ,  i ,  j  ,  k)  /qq  (2 ,  i,  j  ,  k)  )  *57 . 3 
25  continue 
umper=um/62 

write (7, 300) umper , (v(l) ,1=1,5) 
write (8, 300) umper, (ang(l) ,1=1,5) 


30  continue 


write(4,310)  1 , (resd ( 1 , 1 ) , 1=1 , 5) 
do  40  j=2,nres 
it=10* { j -1) 

40  write (4, 310) it, (resd(j ,1) ,1=1,5) 


300  format (6f 8 . 3 ) 

310  format (i5 , 5 (lx, elO . 3) ) 

610  format ( '  *****  warning  *****',/, 

1  '  im,  jm,  )an,  read  from  input', 3i5,'  do  not 
match' , /, 

2  '  im,  jm,  km,  read  from  restart  file',3i5) 
stop 

end 
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This  program  converts  the  x,y  data  point  into  r,z,  theta 
oordinates  for  TCGRID: 

dimension  x ( 5 ) , y ( 5 ) , r ( 5 ) 

data  r/ 3. 941, 4. 033,4.24 0,4. 446, 4. 539/ 

do  10  theta=- .105, .105, .014 

do  20  1=1, 5 

rad=r ( 1 ) 

y  ( 1 )  =  ( rad* sin ( theta) ) *25 . 4 
X ( 1 ) = ( rad* cos { theta) -3.895) *25.4 
20  continue 

write (9, 100) (x(l),y(l).l=l,5) 

10  continue 
100  f ormat { lOf 7 . 3 ) 
stop 


APPENDIX  B.  RVC3D  USER'S  GUIDE 


INTRODUCTION 

RVC3D  is  a  computer  code  for  analysis  of  three- 
dimensional  viscous  flows  in  turbomachinery.  Some  printed 
output  is  available,  but  graphics  software  will  be  needed 
for  examining  the  solution  files.  Restart  files  are  in  the 
standard  qfile  format  for  the  plot3d  or  fast  graphic  codes 
developed  at  NASA  Ames  Research  Center. 

C-type  grids  are  used  to  give  good  resolution  of  blade 
leading -edges  and  wakes.  Grid  input  is  in  standard  plot3d 
xyz-file  format,  so  any  C-grid  generator  can  be  used. 

However  two  grid  codes  have  been  developed  specifically  for 
use  with  rvc3d:  stack  and  tcgrld.  STACK  reads  a  2-D  grid 
generated  by  the  grape  code,  and  generates  a  3D  grid  for  a 
linear  or  annular  blade  row  by  stacking  the  2-D  grid 
spanwise.  tcgrld  is  a  general  3-D  C-  or  H-grid  generator 
for  turbomachinery.  It  reads  annulus  and  blade  geometry  in 
either  merldl  format  or  NASA  Lewis  compressor  design  code 
format.  It  generates  C-type  grids  at  several  spanwise 
locations  using  a  version  of  the  grape  code,  then  reclusters 
the  grids  spanwise. 
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•  •  •  •  •  •  •  • 


STASDENT  SETUP  AND  EXECUTION 


Parameter  statements  are  used  to  make  redimensioning 
simple.  The  parameter  statement  is  located  throughout  the 
code  and  must  be  modified  to  the  dimensions  of  the  input 
grid.  Using  a  vi  editor,  this  can  be  accomplished  by  typing 
the  following  command  at  the  editor  prompt: 

: 1 , $  s/ni=old, nj  =old, nk=old/ni=new, nj=new, nk=new/g 

the  j  component  must  be  one  greater  than  the  j  dimension  of 

the  grid. 

After  the  parameter  statement  is  changed  the  code  must 
be  recompiled,  with  the  following  statement: 

>fc  -03  rvc3d.f  -o  rvc3d 

The  input  file  for  rvc3d  must  also  be  modified  to  set  up 
for  the  flow  conditions.  The  steps  for  this  can  be  found  in 
the  program  documentation. 

Rvc3d  is  run  as  a  standard  unix  process: 

>rvc3d  <  input  >  output 

It  is  recommended  that  this  input  be  put  in  a  .com  file 
and  then  submitted  as  a  batch  job.  For  example: 

>batch  <  rvc3d.com 


APPENDIX  C.  COBRA  PROBE  MEASUREMENTS 


This  table  shows  the  flow  angles  as  measured  by  the 


cobra  probe. 


%  CHORD 

10% 

25% 

50% 

■75% 

90% 

1 

0.0 

62 

65 

62 

62.5 

66 

0.05 

62.5 

65 

63 

63 

65 

1  w 

0.1 

63 

65.5 

63 

63 

65 

0.15 

63 

65 

63 

63 

64.5 

0.2 

62.5 

65 

63 

62.5 

63 

m 

0.25 

61 

63 

62.5 

62 

64.5 

0.3 

61 

62.5 

62 

62 

65.5 

0.35 

61 

64 

62 

62 

73 

0.4 

61 

65 

62 

62 

74 

• 

0.45 

61.5 

65.5 

62 

62 

■73 

0.5 

61.5 

65 

61.5 

61.5 

7 

0.55 

62 

65 

62 

62 

71 

• 

0.6 
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65 

70 

0.65 

63 

64.5 

63 

63 

69 

0.7 

63 

64.5 

63 

63 

68 
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63 
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63 

63 

67 

• 

0.8 

63 

65 

63 

63 

66 

0.85 
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65 

63.5 

63 

65.5 

0.9 
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64.5 
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63.5 

65 

0.95 
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64 
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65 

• 

1.00 
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65 

64 

62.5 

65 

1.05 

63 

65 

63 

63 

66 

1.1 

63 
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63 

63 

65.8 

• 

1.15 

63 

65 

63 

63 

66 

1.2 

63 

65 

63 

63 

65.5 

1.25 

62.5 

63 

62.5 

62.5 

65.5 

1.3 

62 

62.5 

62 

62 

64 

• 

1.15 

1.2 

1.25 


FLOW  ANGLE  AT  25%  SPAN 


%  CHORD 


FLOW  ANGLE  AT  50%  SPAN 


FLOW  ANGLE  AT  15%  SPAN 


Cobra  Probe  Calibranon 


Velocities  on  Selected  Radii 
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%  CHORD 


The  one  component  LDV  survey  data,  for  mid-span  is  tabulated 
below: 


Data  Point 

Velocity  (M/S' 

Turb.  Int.  (%) 

1 

76.848 

6.602 

2 

wmmm 

6.989 

3 

76.067 

7.328 

4 

76.246 

6.872 

5 

76.513 

7.214 

6 

75.873 

7.841 

7 

75.908 

7.601 

8 

76.185 

7.265 

9 

76.423 

7.190  1 

10 

76.256 

7.015 

11 

76.276 

7.877  1 

76.103 

7.766 

13 

75.889 

7.498 

14 

76.416 

7.249 

15 

76.294 

6.992 
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